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Resumen 


Se desarrolla un proceso metodológico para el estudio de la variación espacio- 
temporal de la precipitación pluvial en México, considerando información histórica 
de al menos veinte años a partir de 2 165 estaciones climatológicas pertenecientes al 
Servicio Meteorológico Nacional; se eliminaron datos atípicos y fueron estimados 
con ClimGen, siempre y cuando la información faltante no sobrepasara el 25% 
del total. Se probó bondad de ajuste con respecto a la distribución gamma de dos 
parámetros en los periodos anual y mayo-octubre; este último, por registrar el 80% 
de la precipitación pluvial anual en México. Se observó que en 1 786 estaciones del 
periodo anual y en 1 727 del periodo mayo-octubre se ajustaron a la distribución 
gamma. Se evaluaron las técnicas de interpolación espacial: kriging, cokriging, 
inverso de la distancia y thin plate smoothing spline; este último registró menor 
cuadrado medio del error, y se consideró en la interpolación de los parámetros 
de forma y escala, para generar mapas de precipitación pluvial al 70, 80 y 90% 
de probabilidad. Finalmente, se analizó la frecuencia, duración e intensidad de 
periodos húmedos y secos en México, considerando la clasificación del Índice 
Estandarizado de Precipitación (SPI),, el cual se interpoló para generar mapas 
raster anuales en el periodo 1961-2003, con resolución de 900 metros; se observó 
incremento en la duración de las sequías y disminución en su intensidad. Esta 
investigación pretende coadyuvar a mejorar la modelación del ciclo hidrológico, al 
estudio del cambio climático y a la delimitación de zonas vulnerables en México. 
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Introducción les del ciclo hidrológico (Gourou et al., 1977). 

Sin embargo, la cambiante cantidad, duración, 
La precipitación pluvial abastece de agua dulce frecuencia e intensidad con que ocurre, origina 
al planeta y es una de las variables fundamenta- que sea uno de los fenómenos atmosféricos más 
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difícil de modelar y predecir, además de que 
es uno de los hidrometeoros más afectados 
por eventos naturales, como los llamados 
fenómenos de “El Niño” o de “La Niña”. En 
México, el 80% de la lluvia precipitada en un 
año hidrológico ocurre de mayo a octubre. En 
este periodo se siembran cerca de 12 millones 
de hectáreas de temporal con cultivos cíclicos 
(SIAP, 2008). 
superficie se pierde, entre otras causas, por la 


Normalmente 10% de esta 
falta o el exceso de lluvia. Cabe señalar que este 
ciclo de siembra conocido como de primavera- 
verano de temporal (PV) aporta el 50% del 
volumen total de la producción agrícola. Es 
importante destacar que en México, la toma 
de decisiones respecto a en qué periodos y 
superficies sembrar se sustenta en muchos de 
los casos en el cálculo de la lluvia promedio 
histórica precipitada en un cierto lugar, o bien 
se calcula su probabilidad y anomalía con base 
en la distribución normal, lo cual matemática 
y estadísticamente es posible; sin embargo, 
los resultados generalmente son incorrectos. 
Wilks (1995) menciona que las series de tiempo 
de precipitación se ajustan a una distribución 
de probabilidad tipo gamma. Haddad y 
Rahman (2005) verificaron mediante pruebas 
de Kolgomorov-Smirnov y Anderson-Darling 
que los datos de precipitación registrados en 
el 65% de las estaciones del estado de Victoria, 
Australia, se ajustaron a la distribución gamma. 
García et al., en 1973, lo confirmaron para 
México. No obstante, estos autores utilizaron 
el método de momentos para el cálculo de los 
parámetros de forma y escala, el cual, en sentido 
técnico, puede ser “ineficiente”, ya que no 
hace uso máximo de la información contenida 
en la base de datos. La consecuencia práctica 
de esta “ineficiencia” de acuerdo con Wilks 
(2006) es que los valores pueden ser erráticos 
o “innecesariamente” variar de una muestra 
a Otra. Una mejor aproximación para fijar el 
valor de los parámetros de esta distribución 
es usar el método de máxima verosimilitud, el 
cual requiere un proceso iterativo que sólo es 
práctico mediante el uso de la computadora. 
Cabe señalar que García et al. (1973) no 


presentan en sus resultados los valores críticos 
de alguna prueba estadística de bondad de 
ajuste (Anderson Darling, Kolmogorov y 
Smirnov, o Chi-cuadrada), lo que permitiría 
valorar de una manera estadísticamente acep- 
table el ajuste de los datos de precipitación a 
una distribución gamma. Por otra parte, una 
vez valorada la bondad de ajuste de los datos 
puntuales, es necesario conocer un método de 
interpolación espacial que permita calcular 
valores desconocidos con el mínimo error 
posible. Díaz et al. (2008) evaluaron diferentes 
interpoladores espaciales, concluyendo que el 
que menor cuadrado medio del error presentó 
al interpolar datos de precipitación fue el de 
thin plate smoothing spline, desarrollado por 
Hutchinson (2004) e incluido en el programa 
de cómputo Anusplin. Es importante destacar 
el Índice Estandarizado de Precipitación, 
SPL, por sus siglas en inglés (Standardized 
Precipitation Index), desarrollado por Mckee et 
al. (1993), para determinar si en una región y 
en un periodo determinado hay déficit o exceso 
de precipitación respecto a las condiciones 
normales. Hayes et al. (1999), Tadesse et al. 
(2004) y Komuscu (1999) destacan que el 
método SPI por su rapidez, gran aproximación 
en el análisis de sequía, su simplicidad y 
requerimiento mínimo de datos lo hacen un 
método ideal para utilizarse. Fue diseñado para 
ser un indicador que reconoce la importancia 
de la escala del tiempo en el análisis de la 
disponibilidad y uso del agua. El SPI puede 
ser fácilmente interpretado por un proceso 
simple de medias móviles (Guttman, 1998). 
Es importante señalar que en su algoritmo de 
cálculo se asume que los datos de precipitación 
se ajustan a una distribución gamma. Con base 
en lo anterior, el propósito del presente trabajo es 
aplicar técnicas estadísticas para la evaluación 
de la bondad de ajuste de la precipitación a 
la distribución gamma, evaluar métodos de 
interpolación espacial y aplicar los resultados 
obtenidos para la evaluación de periodos 
de sequía y excesos de humedad, utilizando 
para esto el SPI. Lo anterior permitirá hacer 
cálculos más aproximados de la incidencia 
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probabilística de la precipitación en México, 
generar mapas continuos de probabilidad, 
observar la variación de los parámetros de 
forma y escala de la distribución gamma, así 
como delimitar espacialmente las zonas de 
nuestro país afectadas frecuentemente por 
periodos de sequía o exceso de humedad. Por 
lo anterior, en el presente trabajo se muestra el 
desarrollo de una metodología para el análisis 
espacial y temporal de la precipitación pluvial 
en México, considerando técnicas de simula- 
ción y cómputo iterativo. 


Materiales y métodos 


Los datos de precipitación provienen de 
la Red de Monitoreo Meteorológico de 
la Comisión Nacional del Agua (CNA) y 
fueron consideradas sólo las estaciones que 
contenían al menos veinte años de información 
entre 1961 y 2003, y que la información 
faltante no sobrepasara el 25% del total. 


Prueba de bondad de ajuste y cálculo de 
parámetros de la distribución gamma 


La función probabilística de densidad (fpd) de 
muchas variables atmosféricas son claramente 
asimétricas y sesgadas a la derecha (Wilks, 
1995), algunos ejemplos comunes son la 
cantidad de precipitación pluvial que ocurre 
en un lugar determinado y la velocidad del 
viento. Ajustar estas variables a distribuciones, 
como por ejemplo la Gaussiana, podría llevar a 
resultados incorrectos, al modelar o cuantificar 
cierta probabilidad de ocurrencia. Como se 
ha asentado anteriormente, una distribución 
continua que se ha observado ajusta bien a datos 
de la precipitación pluvial es la distribución 
gamma. Quizá fueron Barger y Thom (1949), 
y Thom (1949), citado por Mosiño y García en 
1981, quienes primero llamaron la atención 
de utilizar la función de distribución gamma 
de dos parámetros para representar datos de 
lluvia. Esta distribución se define por la fpd: 


jay nel) 


O) 


x,a,p>0 


Donde x es un valor de precipitación 
observada, y a y B son los parámetros de forma 
y escala de la distribución. La función T(0.) está 
definida por la integral: 


T(a)= fe"! ar 
0 


Ésta debe ser evaluada numéricamente 
O aproximada utilizando valores tabulares 
provistos en tablas (Wilks, 1995; Sánchez 
et al., 2008). Una de las dificultades de la 
distribución gamma es el cálculo de sus 
parámetros y para ello varios autores proponen 
diferentes formas para estimarlos, uno de 
los mejores es el de máxima verosimilitud, 
el cual requiere un proceso computacional. 
Wilks (1995) describe las aproximaciones 
para este propósito, el primero de ellos es el 
señalado por Thom (1958) y el segundo es una 
aproximación polinomial para el parámetro 
de O propuesto por Greenwood y Durand 
(1960). Gran parte del análisis realizado en 
este trabajo se fundamenta en el trabajo de 
Cheng (2006), quien desarrolló una serie de 
algoritmos y programas en Visual Basic in- 
corporados en hojas de cálculo, con los cuales 
fue posible, mediante el método de máxima 
verosimilitud (MLE), estimar los parámetros 
o. y B de la distribución gamma, sus intervalos 
de confianza, funciones de distribución acu- 
mulada, empírica y de probabilidad, así 
como realizar la prueba de bondad de ajuste 
mediante remuestreo “Bootstrap” (por su 
nombre en inglés) empleando la prueba de 
Anderson-Darling, con su correspondiente 
valor crítico y valor p. 


Evaluación de métodos de interpolación 


Para llevar a cabo investigaciones sobre el 
impacto del clima en diferentes ámbitos 
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es necesario entender la variación de las 
condiciones climáticas (Hartkamp et al. 
1999; Hong et al., 2005). La interpolación de 
variables climáticas ha sido ampliamente 
estudiada en todo el mundo y se han 
empleado varias técnicas para generar 
superficies continuas tanto de precipitación 
como de temperatura (Dressler et al., 2000). 
El proceso de interpolación se considera de 
suma importancia, en virtud de que, en un 
plano, se asigna un valor desconocido a un 
punto georreferenciado, partiendo de valores 
conocidos cercanos a ese punto. Algunos 
esfuerzos reportados son los de Hartkamp 
et al. (1999), Téllez et al. (2003), Jiménez et 
al. (2004) y Treviño et al. (2002), donde el 
común denominador de estos estudios fue la 
aplicación de técnicas de interpolación para la 
obtención de superficies continuas de distintos 
temas climáticos. Por otro lado, Descroix et al. 
(1997) realizaron un estudio en la cuenca alta 
del río Nazas, donde se encontró correlación 
de la precipitación con la altitud (0.85) y 
distancia al océano (-0.85), indicando que 
la precipitación se incrementa con la altura 
y decrece a medida que se aleja del océano, 
explicando un 85% la variación de la cantidad 
de lluvia. De estos trabajos se puede concluir 
que el método de interpolación empleado 
debe ser en función de las características de la 
variable de interés, su continuidad espacial, la 
superficie geográfica sobre la cual se va aplicar 
y el número de datos disponible. 


Resultados 


De las 5 088 estaciones climáticas analizadas, 
sólo 2 165 cumplieron los criterios de cali- 
dad y cantidad de datos establecidos con 
anterioridad. Dichas estaciones se integraron 
a un Sistema de Información Climática (Díaz y 
Sánchez, 2007) para facilitar su procesamiento 
y análisis. Los datos atípicos más comunes 
detectados con R-ClimDex (Zhang y Yang, 2004) 
fueron temperaturas máximas por debajo de 
las mínimas, temperaturas mínimas rebasando 
las máximas y precipitaciones negativas. En 
caso necesario, los datos faltantes, como ya 
se indicó, fueron estimados con el ClimGen 
(Nelson, 2003). En cuanto a la evaluación de 
la bondad de ajuste, las rutinas desarrolladas 
por Cheng (2006) resultaron ser prácticas y 
rápidas para calcular los parámetros de forma 
(a) y escala (PB), valores críticos, valores-p, e 
intervalos de confianza de dichos parámetros. 
La aplicación iterativa de estas rutinas 
permitió estimar los parámetros y estadísticos 
de la prueba de bondad de ajuste para los 
periodos mayo-octubre y anual. El cuadro 1 
corresponde a los resultados del análisis de 
bondad de ajuste a los datos de precipitación 
anual; sólo se listan algunas estaciones con 
los parámetros y estadísticos de prueba 
resultantes, y se encuentran ordenadas por 
el valor-p. En la parte superior se encuentran 
las estaciones con probabilidades altas de que 
los datos de precipitación se ajusten a una 


2 

SS] 
3 Cuadro 1. Parámetros y estadísticos de prueba de la bondad de ajuste a la distribución gamma para el periodo anual. 
E 

2 
z Parámetro Al 

cti 
E Estación Años | Pp-Mín | Pp-Máx | Pp-Med | Forma | Escala | AD pa 12 | valor Pp 
3 (a) (B) 

e Bustamante 34 279 864 513 13.98 36.65 | 0.14 0.76 0.98 
S 

al Yurécuaro 85 486 Ala 763 22.61 33.74 | 0.15 0.75 0.97 
? Sto. Domingo Teojomulco 36 1229 3152 2 095 33.67 62.22 | 0.69 0.74 0.07 
» Chicontepec Tejada 40 1216 4 163 1947 11.30 727379 |00:98 0.74 0.02 
$ San Miguel del Monte 23 183 1949 1018 4.44 229.03 | 1.00 0.74 0.01 


« vj30/0u32], 


Díaz-Padilla et al., Variación espacio-temporal de la precipitación pluvial en México: una aproximación a la evaluación de impactos 


distribución gamma y en la parte inferior 
se tienen las que registraron menor valor-p 
(menor a 0.05), lo que conduce a rechazar la 
hipótesis de que los datos siguen el patrón de 
la distribución antes mencionada. 

Posterior a la evaluación de la bondad 
de ajuste se interpolaron los parámetros de 
forma (0) y escala (PB), mediante el programa 
Anusplin, que con base en el trabajo de Díaz et 
al. (2008) resultó con menor cuadrado medio 
del error de predicción, al compararle con los 
interpoladores espaciales: kriging y cokriging. 
La figura 1 muestra el comportamiento 
geográfico del parámetro Q;; la figura 2 muestra 
el comportamiento del parámetro PB, cuya 
información proviene de la precipitación 
pluvial anual registrada en 1 786 estaciones. 

En cuanto a la precipitación anual acu- 
mulada y mayo-octubre, el parámetro 0, 
representado en la figura 1, tiende a formar 
un patrón simétrico en la región de la Sierra 
Madre Occidental y no en la región de la Sierra 


Madre Oriental ni en la región centro sur y sur- 
sureste del país; en tanto que en regiones como 
la península de Baja California, Sonora, y gran 
parte de la región noreste y norte-centro, tiende 
a presentar fuertes sesgos. 


Monitoreo de sequías y excesos de 
humedad en México mediante el SPI 


La demostración de algunas de las aplicaciones 
basadas en los resultados obtenidos y herra- 
mientas desarrolladas en el presente trabajo 
se describen a continuación. El cálculo del SPI 
sirve para determinar si en una región y en 
un periodo determinado hay déficit o exceso 
de precipitación en cuanto a las condiciones 
normales de ocurrencia. Fue diseñado para 
ser un indicador que reconoce la importancia 
de la escala del tiempo en el análisis de la 
disponibilidad y uso del agua. El cálculo del 
SPI se basa en el ajuste de la serie histórica de 
precipitación a la distribución gamma que se 
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Figura 1. Delimitación geográfica del parámetro a de la distribución de probabilidad gamma en el territorio nacional, 


correspondiente a la precipitación pluvial anual registrada en 1 786 estaciones. 
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Figura 2. Delimitación geográfica del parámetro f de la distribución de probabilidad gamma en el 


territorio nacional, correspondiente a la precipitación pluvial anual registrada en 1 786 estaciones. 


transforma a una distribución de probabilidad 
normal estandarizada. De modo que el valor 
medio del SPI es cero y varianza uno, lo que 
indica que los valores positivos son mayores a 
la mediana y los negativos son menores a ésta 
(Crespo, 2005). 

Para ello se estiman los parámetros a 
y B de la distribución para cada escala de 
tiempo de interés. Estos parámetros se 
utilizan para encontrar las probabilidades 
acumuladas de un evento de precipitación en 
un periodo determinado de acuerdo con los 
datos disponibles de cada estación. Dichas 
probabilidades acumuladas se transforman en 
una variable estándar con media igual a cero y 
varianza igual a 1, lo cual es el valor del índice 
SPL donde tiene la ventaja de ser calculado 
para distintos periodos de tiempo, por ejemplo 
1, 3, 6, 12 o 72 meses. Una descripción paso a 
paso de su estimación puede consultarse en el 


trabajo realizado por Giddings y Soto (2005). 


A continuación, el cuadro 2 presenta la ca- 
tegorización de los valores de SPI de acuerdo 
con los valores obtenidos después del proceso, 
dicha clasificación es la que utiliza el National 
Drought Monitoring Center, y algunas otras 
instituciones alrededor del mundo. 

Se realizó el estudio para la estación 
climatológica de Pánuco, Veracruz, México, 
donde se interpretan sus valores de SPI 
para las escalas de tiempo de 6 y 12 meses. 
En lo referente a la escala de tiempo de seis 
meses, la intensidad de eventos de humedad 
se presentó a principios de 1968, poco antes 
de medio año de 1982, y principios de 1992; 
este último es uno de los periodos con mayor 
duración, puesto que abarca de finales de 
1990 hasta mediados de 1992; otro evento 
importante fue el de 1973 hasta principios 
de 1975; en ambos periodos se nota un ligero 
descenso, para posteriormente incrementarse 
nuevamente a muy húmedo y extremadamen- 
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Cuadro 2. Clasificación de los valores del SPI de acuerdo 
con el National Drought Monitoring Center. 


Valor de SPI Categorías de sequía 
20 mayor Extremadamente húmedo 
15a1.99 Muy húmedo 
1.0a 1.49 Moderadamente húmedo 
-0.99 a 0.99 Cercano a lo normal 
-1.0 a -1.49 Moderadamente seco 
Iacio9 Muy seco 
-2 0 menor Extremadamente seco 


te húmedo. Respecto a las condiciones de 
sequía, se distinguen cinco eventos de cate- 
goría extremadamente seco, sin embargo, los 
años en que se hace más notable son 1971, 1980 
y 1982. En cuanto a intensidad y duración, se 
destacó el evento se sequía ocurrido a finales 
de 1962 hasta principios de 1964. En cuanto 
al periodo anual, se observó que el evento 


más húmedo se presentó a principios de 1992, 
con una duración desde finales de 1990 hasta 
finales de 1992, presentando dos “picos” 
extremadamente húmedos; en la categoría 
de muy húmedo es donde se hace presente 
la mayor duración, ya que abarca a partir de 
medio año de 1973 hasta 1977. En relación con 
la sequía, se apreció un evento extremadamen- 
te seco a mediados de 1965 y otro con muy 
poca presencia en 1983; en intensidad y 
duración de algún evento de sequía ocurrió 
a finales de 1961 hasta finales de 1965, donde 
alcanzó la categoría de extremadamente 
seco. En el lapso de 1994 a 2001 se presentó 
un evento clasificado como moderadamente 
seco, también considerado como de mayor 
duración. Por último y en términos generales, 
es importante hacer mención que la frecuencia 
de la sequía en los periodos de 6 y 12 meses se 
hace cada vez más presente en sus diferentes 
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Figura 3. Valores del SPI generados para diferentes escalas de tiempo en la estación climatológica Pánuco del 


estado de Veracruz, para el registro histórico de la precipitación del periodo 1961-2003. 
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categorías, haciéndose más ausentes los 
periodos de humedad (figura 3). 

Otra de las ventajas del cálculo de los 
valores del SPI es que permite la comparación 
de escalas temporales de interés particular, 
por ejemplo, en la figura 4 se presentan los 
valores de SPI para el periodo mayo-octubre 
de Pánuco, Veracruz, para el periodo 1961- 
2003. 

De la figura 4 se puede apreciar la frecuen- 
cia de eventos extremos manifestados en dos 
de cada 27 años; las condiciones del temporal 
van a ser de sequía extrema y de condiciones 
muy húmedas en la misma proporción. Res- 
pecto a la asociación de sequías extremas 
con años Niño o Niña, se puede indicar que 
1982 fue uno de los años más desastrosos de 
la agricultura en México (Magaña y Morales, 
2004) y, como se aprecia en la gráfica, el valor 
del SPI sí captó la presencia del fenómeno 
en 1982; sin embargo, la sequía de 1994 no 
se reporta como evento Niño o Niña, por lo 
que se requerirían estudios más profundos al 
respecto. Magaña y Morales (2004) señalan 
que el evento El Niño más intenso fue el 
ocurrido entre los años 1991 y 1995, el cual 
fue registrado perfectamente por el SPI como 


sequía extrema. 


Incertidumbre climática y política 
pública 


En la definición e implementación de políticas 
públicas dirigidas a promover el desarrollo 
rural sustentable y mejorar las condiciones 
de vida de los productores del campo 
mexicano es necesario, como ha quedado de 
manifiesto en el presente trabajo, cuantificar 
y delimitar el riesgo climático. Esto significa 
considerar los diagnósticos y pronósticos cli- 
máticos realizados psara diferentes variables 
atmosféricas que impactan en el desarrollo y la 
productividad de los cultivos, la conservación 
de los recursos naturales y el bienestar social 
en general. Este esquema de análisis ayudaría 
a evitar impactos negativos significativos que 
redundarían en elevados costos económicos y 
sociales para el país. En la figura 5 se muestra 
la distribución espacial del SPI durante el 
año 1982, notándose que más del 50% de la 
superficie nacional se vio afectada por sequía 
severa, demeritando las bondades hipotéticas 
de los programas nacionales promulgados en 
esa época. Las causas fueron desastrosas para 
el país, en donde cerca del 40% de la agricultura 
de temporal sembrada con maíz fue declarada 
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Figura 4. Frecuencia, intensidad y severidad de eventos de sequía y humedad para la estación de Pánuco, Veracruz. 


Díaz-Padilla et al., Variación espacio-temporal de la precipitación pluvial en México: una aproximación a la evaluación de impactos 


-115 -110 -105 


+ + 
Índice Estandarizado de Precipitación 
Periodo de tiempo, 12 meses (enero-diciembre), 
correspondiente al año 1982 


+ 


600 Kilómetros 


-115 -110 -105 


-100 -95 -90 


_ m + 
Instituto Nacional de Investigaciones 
Forestales, Agrícolas y Pecuarias 


Categoría de la sequía sPI e 
+2.01 a +3.00 | 0.00 
+151a+2.00| 0.00 
+101a+151| 071 
+051a+1.00| 7.54 
-0.49 a +0.50 | 45.79 
0.99 2-0.50 | 23.45 
-149a-1.00 | 17.97 
-199a-150 | 4.54 
-3.00a-2.00 | 0.00 


E533) Extremadamente húmedo 
EZ Muy húmedo 

E Húmedo 
ES5Ligeramente húmedo 
Normal 

YY Ligeramente seco 
E Seco 

EXA Muy seco 

55] Extremadamente seco 


-100 -95 -90 


Figura 5. Distribución espacial Índice Estandarización de Precipitación, donde se ha reportado condiciones de El Niño. 


como siniestrada a causa, principalmente, de la 
intensa sequía. 

Con la idea de concentrar los resultados 
obtenidos se presenta el cuadro 3, que mues- 
tra un concentrado de cómo han variado 
través del tiempo las condiciones de sequía y 
humedad en el territorio nacional. 

De esta manera puede apreciarse que para 
1961, cerca del 73% del territorio nacional se 
ubicó en la categoría de condiciones normales 
y cerca del 17% del mismo en condiciones de 
ligera sequía. 

La información de la figura 6 proporciona 
un reflejo del comportamiento histórico de 
los eventos de sequía, los cuales han ido en 
aumento a partir de la década de los noventa, 
en tanto que las condiciones de humedad han 
venido disminuyendo notablemente a partir 
de la misma década. 


Sin duda, estos indicadores pueden con- 
siderarse para contar con apreciaciones 
objetivas de las probabilidades de éxito 
en acciones con impacto social y adecuar 
la toma de decisiones en consecuencia, al 
considerar los distintos escenarios climáticos 
que pueden ocurrir en el sector agropecuario 
nacional. 


Discusión 


Para la integración de series climáticas de 
calidad y sin datos faltantes, se sugiere el uso 
del software R-ClimDex como una alternativa 
para valorar la calidad y el volumen de datos 
de las series históricas diarias de precipitación. 
ClimGen mostró ser un generador climático 
confiable para la estimación de datos faltantes 
diarios de precipitación, ya que respeta los 


Tecnología y 


Ciencias del Agua, vol. UL. núm. 4, octubre-diciembre de 2011 


Díaz-Padilla et al., Variación espacio-temporal de la precipitación pluvial en México: una aproximación a la evaluación de impactos 


Cuadro 3. Porcentaje de la superficie incidente de la república mexicana acorde con las diferentes clasificaciones 
del Índice Estandarizado de Precipitación (SPI) anual, para el periodo 1961-2003. 


Extremadamente Muy Ligeramente Ligeramente Muy Extremadamente 
Año Húmedo Normal Seco 
húmedo húmedo húmedo seco seco seco 
1961 0.00 0.00 0.34 8.20 72.93 16.79 1.61 0.12 0.00 
1962 0.00 0.00 0.61 2.74 47.39 32.87 16,23 0.15 0.00 
1963 0.00 0.00 L97 10.60 64.44 19.98 2.93 0.08 0.00 
1964 0.00 0.00 0.01 6.03 62.61 21.56 8.85 0.95 0.00 
1965 0.00 0.19 1.14 8.75 64.35 21.85 3,72 0.00 0.00 
1966 0.00 0.00 4.77 24.99 66.23 2.50 0.81 0.41 0.29 
1967 0.35 4.60 14.91 22.66 44.60 11,89 1.51 0.00 0.00 
47.50 8.76 4.36 220 0.34 
40.78 21,52 12.76 3.21 0.07 
62.68 20.53 5.26 0,39 0.15 
51.88 18.59 4.77 0.61 0.00 
62.54 11.48 0.40 0.00 0.00 
45.23 14.23 0.19 0.00 0.00 
61.15 16.66 aa) 1.25 0.17 
56.96 2712 9.24 0.08 0.00 
47.39 4.60 1.30 0.00 0.00 
B5.35 1410 9.76 0.59 0.00 
59.47 2.03 0.00 0.00 0.00 
47.96 27.30 1.97 1.88 0.00 
64.04 17,37 6.03 0.33 0.00 
35.18 1.63 0.00 0.00 0.00 
45.79 23.45 197 4.54 0.00 
56.30 11.45 1.60 0.00 0.00 
40.33 2.44 0.91 0.00 0.00 
76.87 10.14 0.01 0.00 0.00 
48.17 NS 2.61 0.01 0.00 
39.79 31.88 8.93 0.18 0.00 
71.19 13.00 Ur 0.84 0.04 
53.34 33.36 9.01 La 0.07 
44.57 12.30 2.80 0.18 0.00 
E 50.44 939 2.23 0.04 0.00 
3 47.04 1.34 0.00 0.00 0.00 
bs 66.23 9.92 0.63 0.00 0.00 
E 44.09 29.04 10.68 | 573 0.00 
3 41.81 29.63 12.04 0.00 0.00 
É 62.27 22.74 8.60 0.00 0.00 
3 78.42 10.00 0.14 0.00 0.00 
- 48.08 31.49 10.80 0.12 0.00 
E 39.44 28.02 14.71 4.47 0.54 
= 64.94 21.56 4.77 0.60 0.00 
_ 60.77 30.49 6.68 0.15 0.00 
El 63.03 20.55 3.90 0.00 0.00 
$ 2003 0.00 0.00 2.34 11.74 63.14 17.96 4.75 0.07 0.00 
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Figura 6. Comportamiento cronológico de las superficies registradas bajo condiciones de humedad y sequía en México. 


datos existentes y sólo genera los faltantes; 
además, resultó más preciso al compararlo 
mediante pruebas estadísticas, con otros 
generadores climáticos. Respecto a la bondad 
de ajuste de la precipitación a la distribución 
gamma, la distribución de probabilidad teórica 
gamma de dos parámetros es la adecuada para 
la modelación probabilística de la precipitación 
pluvial en México tanto para los registros 
anuales como para los del periodo mayo- 
octubre, y se propone su uso en los cálculos 
de anomalías climáticas y la predicción de 
eventos extremos. No se observó algún patrón 
geográfico, climático o espacial que influyera 
en la bondad de ajuste de esta distribución, 
por lo que podrá utilizarse en cualquier clima 
o espacio geográfico de México. El método 
adecuado de interpolación espacial de la 
precipitación pluvial es el thin plate smoothing 
spline, dado que registró el menor cuadrado 
medio del error de predicción al compararlo 
con algunos otros métodos de interpolación 
(kriging, inverso de la distancia al cuadrado 
-IDW- y cokriging). Se recomienda usar 


los mapas interpolados de los parámetros 
de forma y escala de la distribución gamma 
para modelar el ciclo hidrológico y el cambio 
climático, y para el cálculo probabilístico de 
eventos extremos en México. El SPI resultó 
ser una herramienta sencilla de calcular e 
interpretar, y captó en forma eficiente los 
diversos periodos de humedad y sequía 
en México. Se deduce que con las series 
históricas calculadas e interpoladas del SPI es 
posible monitorear puntual y espacialmente 
la intensidad, frecuencia y duración de las 
sequías o exceso de humedad en México. Al 
concentrar los resultados de la interpolación 
de las series históricas del SPI, se encontró 
que el periodo 1993-2003 se caracteriza por 
un incremento en la duración de las sequías y 
una baja en la intensidad de las mismas. Al no 
considerar los avances de los conocimientos 
científicos y no aplicarlos a la definición de 
políticas públicas, los tomadores de decisiones 
podrían incurrir en errores considerables, 
tal es el caso del lanzamiento del programa 
Sistema Alimentario Mexicano (SAM) en 
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1982, para promover la siembra de maíz en 
México, sin considerar que las condiciones 
climáticas del país estaban afectadas por 
la influencia de un año Niño, por lo que el 
porcentaje de siniestro fue uno de los más 
altos de la historia. La delimitación espacial 
del impacto de eventos climáticos como de El 
Niño y La Niña pueden monitorearse a través 
de los mapas interpolados del SPI y deben 
ser la base para el estudio de la variabilidad, 
cambio climático y zonas de vulnerabilidad en 
México. Finalmente, algunos otros resultados 
y conclusiones relevantes pueden consultarse 
con más detalle en el trabajo de Díaz (2009). 


Conclusiones 


Se desarrolló un proceso metodológico pa- 
ra la delimitación espacio-temporal de la 
variabilidad de la precipitación pluvial en 
México. Con base en el ajuste de las series 
históricas de precipitación a la función de 
densidad gamma de dos parámetros, se 
procedió al cálculo del Índice Estandarizado 
de Precipitación, el cual resultó ser un método 
sencillo, preciso y práctico para detectar 
la intensidad, frecuencia y duración de las 
anomalías o variabilidad puntual y espacial 
de la precipitación pluvial en México. Es 
importante destacar que mediante el cálculo 
e interpolación de series históricas del SPI 
se detectó una tendencia al incremento y 
una baja en la intensidad de los eventos de 
sequía en México. Los resultados obtenidos 
en la presente investigación coadyuvarán, 
entre otros aspectos, al modelaje del ciclo 
hidrológico, al estudio del cambio climático 
y a la delimitación de zonas vulnerables en 
México. 
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Abstract 


DÍAZ-PADILLA, G., SÁNCHEZ-COHEN, I., QUIROZ, R., GARATUZA-PAYÁN, J., WATTS- 
THORP, C., CRUZ-MEDINA, LR. € GUAJARDO-PANES, R.A. Space-time variation of rainfall 
in Mexico: an approach for assessing impacts. Water Technology and Sciences (in Spanish). Vol. Il, 
No. 4, October-December, 2011, pp. 51-64. 


This work discusses the development of a process to approach the study of spatial and temporal rainfall 
patterns in Mexico, considering historical information for at least 20 years from 2 165 weather stations 
belonging to the National Weather Service. Outlier data were eliminated and estimates were made 
using ClimGen, provided that the missing information did not exceed 25% of the total. Goodness-of-fit 
was tested with respect to the two-parameter gamma distribution for the annual period and the period 
from May to October, the latter of which registers 80% of the annual rainfall in Mexico. It was noted 
that 1 786 of the stations for the annual period and 1 727 stations for May to October adjusted to the 
gamma distribution. The following spatial interpolation techniques were assessed: kriging, co-kriging, 
inverse distance and thin-plate smoothing spline. The latter showed lower mean square error and was 
considered in the interpolation of the shape and scale parameters to generate precipitation maps for 70, 
80 and 90% rainfall probability. Finally, we analyzed the frequency, duration and intensity of wet and 
dry periods in Mexico considering the classification by the Standardized Precipitation Index (SPD, 
which is interpolated to generate annual raster maps for the period 1961-2003 with a resolution of 900 
meters. An increase was observed in the duration of droughts and a decrease in their intensity. This 
research aims to help improve the modeling of the hydrological cycle, the study of climate change and 
the delineation of vulnerable areas in Mexico. 


Keywords: precipitation, spatial interpolation, goodness of fit, SPI. 
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